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Abstract 



In this note we discuss a supersymmetric (SUSY) D^ x Z^ model which leads to vanishing 
reactor mixing angle ^13 = and maximal atmospheric mixing ^23 = 7''/4 in the lepton sector 
at leading order (LO), due to the preservation of non-trivial distinct D^ subgroups in the 
charged lepton and neutrino sectors, respectively. The solar mixing angle 612 remains unde- 
termined and is expected to be of order one. Since right-handed charged leptons transform 
QQ \ as singlets under 1)4, the charged lepton mass hierarchy can be naturally accounted for. The 

00 ' model predicts inverted mass hierarchy for neutrinos. Additionally, we show that, unlike in 

most of the other models of this type, all vacuum expectation values (VEVs) of gauge sin- 
[■^^ . glets (flavons) can be fixed through mass parameters of the superpotential. Next-to-leading 

^^ ' order (NLO) corrections to lepton masses and mixings are calculated and shown to be under 

^^ . control, especially the corrections to 623 = 7r/4 and ^13 = are of the order of the generic 

expansion parameter e ^ 0.04 and arise dominantly from the charged lepton sector. 
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1 Introduction 

The observed mixing pattern in the lepton sector [T], 

sin^^ia = 0.318lj];[!^^ , sin^ ^23 = 0.50l[]:J? and sin^ ^13 = 0.013+[!:[!^^ (2ct) , (1) 

is weU compatible with special mixing patterns in which ^13 vanishes and 623 is maximal, 

sin 623 = — and sin0i3 = . (2) 

Eq. ([2]) is derived from a /i — r symmetric neutrino mass matrix [2] , in the charged lepton mass 
basis. At the same time the solar mixing angle 9i2 is not predicted. It is tempting to assume a 
symmetry to be the origin of such a mixing pattern. 

It has been shown [3] that in case of a non-trivial breaking of a dihedral flavor symmetry in 
the lepton (quark) sector one element Uai3 of the lepton (quark) mixing matrix is given in terms 
of group theoretical quantities only |^ 



\U, 



a/3 1 



TT{ki -k2)i 

COS ' 



n 



(3) 



where n refers to the group index of the dihedral group Dn , j to the index of the two-dimensional 
irreducible representation 2; under which two of the three left-handed fields transform and ki^2 

are associated to the generating elements BA ^ and BA ^ of the Z2 subgroups preserved in the 
charged lepton (down quark) and neutrino (up quark) sectors, respectively. (A and B are the 
generators of the original dihedral group Dn-) Note that ki and ^2 have to be distinct to get a 
non-trivial value for |[/a/3|. Considering the lepton sector, one sees that choosing n = 4 makes it 
possible to achieve e.g. \Ufj,3\ = -j=. Using a specific set of flavons to break L>4 to the subgroup 

Z2 generated by BA ^ in the charged lepton sector, also |C/r3| = -j^ and f/gs = can be enforced, 
thus leading to /i — r symmetric lepton mixings. 

In the model presented here D4 is accompanied by the cyclic symmetry Z5 . Both symmetries 
are only spontaneously broken by fiavon VEVs. We use as framework the Minimal SUSY 
Standard Model (MSSM). It is not the first time that jU — r symmetry is deduced from the group 
Da jUEHnilTO] . However, in these models producing the mass hierarchy among charged leptons 
is usually non-trivial. In contrast, the observed charged lepton mass hierarchy 

nie '■ m^ : rrir ~ e : e : 1 with e ~ A ~ 0.04 , (4) 



with A being the Wolfenstein parameter [TT], is a natural outcome of our model. In order to 
achieve this, it is essential that the right-handed charged lepton fields transform as (different) 
singlets under the fiavor symmetry D4 , instead of being in the representations 1 + 2. Notice 
that also the introduction of a U{1) symmetry is not necessary for ending up with Eq.(j3|), see 
as weh [IMlj- 

At LO, our model leads to ^13 = and ^23 = '?i"/4 with no constraints on the solar mixing 
angle 612, apart from predicting its value generically to be of order one. In the limit of preserved 
subgroups rUfj^ and rrir are generated of the correct order of magnitude. Since rrir arises from a 
non-renormalizable operator small and moderate values of tan /3 are preferred. Neutrinos exhibit 



^For specific realizations see [JHS] 



inverted mass hierarchy and the hghtest neutrino mass m^ fulfills m^ > 0.015 eV. We find at LO 
a strong correlation among 771,3 and the measure of neutrinoless double beta decay \mee\, a rather 
weak one among m^ and tan 9i2 as well as a restricted range for Majorana phases. These results 
are very similar to those found in [4j. We study the model to NLO and show that corrections 
to the predictions ^13 = and ^23 being maximal are of order e and dominated by the charged 
lepton sector. Also 612 undergoes corrections; but these are not particularly interesting, because 
its value is not a prediction of our model anyway. The mass of the electron is generated by 
NLO corrections. The masses of the other charged leptons and neutrinos are slightly corrected 
as well. 

Furthermore, an appropriate construction of the fiavon superpotential allows us to fix all 
flavon VEVs through couplings with positive mass dimension. By choosing the latter to be of 
order eA, where A is the generic cutoff scale of our theory, all flavon VEVs are of order eA. Fixing 
all flavon VEVs by mass parameters is usually not the case in models of such type which thus 
encounter free parameters among the flavon VEVs (equivalent to flat directions in the flavon 
(super)potential), see e.g. [H[5l|T2l[15] . 

The paper is structured as follows: in Section [2] we briefly introduce the group D4 and 
discuss its subgroups. The model is outlined and the LO results for lepton masses and mixings 
are presented in Section [3l Section 2] contains the discussion of the flavon superpotential at LO 
and NLO. In Section [5] NLO corrections to lepton masses and mixings are shown to be well 
under control. We conclude in Section [6l 

2 D4 Group Theory 

The dihedral group D4 has eight (distinct) elements and five irreducible representations denoted 
here as ]^j, i = 1, ..., 4 and 2. All representations are real and only 2 is faithful. L'4 is generated 
by A and B which can be chosen as [16] 



for 2. Note that we have chosen A to be a complex matrix, although 2 is a real representation. 
For (oi, 02)* ~ 2 we thus find that (02, a*)* transforms as 2 under 1^4. The generators A and B 
of the one-dimensional representations are 

11 : A = l, B = l, (6) 

12 : A = l, B = -l, (7) 

13 : A=-l, B = l, (8) 

14 : A=-l, B = -l. (9) 

The character table can be found in e.g. [Ij. A and B fulfill the relations 

A^ = 1 , B^ = 1 and ABA = B . (10) 

The Kronecker products involving Ij are the following 
li xli = ll , ll xli = li for i = 1, ..., 4 , I2 X I3 =14,12 X I4 = I3 and I3 x I4 = I2 . 



For Si ~ Ij and (ai, 02)* ~ 2 we find 



s,a, \ 2 , "^"1 ]-2 , ''"' U 2 and 

S1O2 / V -•5202 / V ''^Sfll / V -■5401 



The four one-dimensional representations contained in 2 x 2 read for (oi, 02)*, (61, 62)* ~ 2 
0162 + ^2^1 ~ il , 01^2 - ^2^1 ~ I2 > o-ih + ^2^2 ~ I3 and aibi - 0262 ~ 14 . 

These formulae are special cases of the expressions, given for dihedral symmetries -D„ with 
general index n, which can be found e.g. in [3lll7]. 

In order to understand how maximal atmospheric mixing and vanishing ^13 arise in our 
model, it is relevant to discuss the subgroups of D4. All its subgroups are abelian: Z2 = Di, Z4 
and D2 = Z2 X Z2- In the following we are only interested in Z2 subgroups generated by BA 
with k = 0, ...,3, because these are preserved in the charged lepton and neutrino sectors at LO. 
Noting 

(BA'=)2 = BA^BA'^ = BA'^-1BA'=-i = ••• = B2 = 1 

shows that BA indeed generates a Z2 symmetry. Since k can take integer values between 
and 3, we find four possible Z2 subgroups of this type. Apart from VEVs of fields transforming 
trivially under D4, a Z2 group given through BA is left unbroken by a non- vanishing VEV of 
a singlet transforming as I3 if fc = 0, 2 holds, and of one singlet transforming as I4 for A; = 1, 3. 
Additionally, it is left intact by fields '(/'i,2 forming a doublet, if their VEVs have the following 
structure 

^^i> ^o.( '~^] . (11) 



(V'2> ; \^ 1 ^ 

As mentioned in the Introduction, in order to get non-trivial mixing the Z2 subgroups preserved 
in the charged lepton and the neutrino sectors have to have different indices ki and k,^. To 
achieve maximal atmospheric mixing, we need e.g. /c/ to be odd, whereas fejy has to be even. 
This is analogous to the constraints found for the indices ki^^, in [Ij. In the following section we 
show that it is phenomenologically irrelevant whether fc; is 1 or 3 as well as whether ki, is or 
2, i.e. the result for the lepton mixing angles is in all cases 6*13 = and 6*23 = 7r/4. In Section 
m a simple fiavon superpotential is constructed which naturally leads to an odd index ki and an 
even index k^. 

3 Outline of the Model and Results at LO 

In this section the transformation properties of all fields (apart from the driving fields responsible 
for the vacuum alignment) under the fiavor symmetry 1)4 x Z5 are presented. Similar to the 
model discussed in [3] the framework is the MSSM and the left-handed lepton doublets Lj 
transform as Ij^ + 2 under D4. In order to accommodate the charged lepton mass hierarchy, 
we assign the right-handed charged leptons e? to the three singlets I2 + I3 + I4 and not to 
ll + 2 as done in [4J . In contrast to the model given in 2] we assign the Z5 charges to the 
matter superfields and hu^d such that hu^d are uncharged under D^ x Z5 and, more importantly, 
we allow the generations of one type, Li and ef , to transform differently under Z^. In this way 
the additional cyclic symmetry Z5 does not only play the role of a symmetry which separates 
the charged lepton and neutrino sectors (at LO), but also plays the role of a Froggatt-Nielsen 
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Table 1: Particle content of the model. Li denotes the three left-handed lepton 
SU{2)l doublets, el the right-handed charged leptons and hu,d the MSSM Higgs 
doublets. The flavons V'e;i,2 and V'i';i,2 only transform under D4, x Z5. Also the 
Froggatt-Nielsen type fields r]\ and rys are gauge singlets and only carry a non-zero 

27ri 

Z^ charge, uj — e ^ is the generating element of Z^. 



symmetry [18]. As one can see, we use for the five different multiplets Li, L^ = (^2,^3), ef 
all five different possible Z5 charges. As in [ij the light neutrino masses arise from the effective 
operator LjLj/i^/A. The flavons relevant for achieving n — t symmetry in the lepton sector at 
LO, are f/'e and ipu which form doublets under Di^. The VEV of ijje is aligned through the flavon 
superpotential in such a way that either the Z2 subgroup generated by BA or by BA'^ is preserved 
in the charged lepton sector, whereas the special form of {'ipu;i,2) preserves either a Z2 subgroup 
arising from B or from B A^ . As we will show it is not relevant for phenomenology which of the 
two subgroups in each sector is chosen. The relevant aspect is the fact that different subgroups 
are preserved. The two additional flavons r/i and % are Froggatt-Nielsen type fields. Since they 
do not transform under D4^ they do not play a role in the preservation of different subgroups of 
1)4. They allow for non-zero (11) and (23), (32) entries in the neutrino mass matrix at LO and 
are relevant for generating the appropriate hierarchy among the charged lepton masses, because 
rrir arises dominantly from L^e^hdipe/^, while the mass of the muon is generated through 
the operator LDc'^hdipeVi/^'^- The mass of the electron originates mainly from one subleading 
operator involving three flavons, Lie^hd'tpe'^uVi / ^"^ ■ Since the tau lepton mass stems from a non- 
renormalizable operator, Vd = (hd) is expected to be of the order of the electroweak scale and 
thus small and moderate values of tan/3 = {hy) /{hd) = Vu/vd are preferred in this model. The 
right order of magnitude of mr can be achieved for ('i/;e)/A ~ e ~ 0.04. The ratio rn,^ : nir ~ e : 1 
requires that also {r]i)/A ~ e. As derived from the flavon superpotential all flavon VEVs can 
be of order eA and thus the electron mass rrie is expected to fulflll rUe : rur ~ e^ : 1 due to 
its origin from a three flavon operator. Apart from generating the electron mass, operators 
involving several flavons (as well as shifts in the flavon VEVs) lead to small deviations from 
the LO result that 623 is maximal and ^13 vanishes. As we show below, these deviations are 
dominated by corrections associated to the charged lepton sector, while all subleading effects in 
the neutrino sector are of relative order e^. All fields appearing in the Yukawa couplings and 
their transformation properties under D4 x Z5 are captured in Table [TJ 

In order to elucidate how maximal atmospheric mixing and vanishing reactor mixing angle 
arise, we first take into account only operators which are suppressed by at maximum 1/A^ and 
in which the fiavon doublet ipi, only couples to neutrinos, whereas V'e couples to charged leptons. 
In this way the Z2 subgroups in charged lepton and neutrino sectors remain exactly preserved. 
We show subsequently that the inclusion of further operators, compatible with the symmetries 
of the model, and shifts in the flavon VEVs only slightly correct the results achieved at this 



level. With the above restrictions the allowed operators in the superpotential aren 

6 C 

Wi= ^ (i2e3/idVe;l " L^elhdi^e,2) + ^ (LaC^/irfVe;! + L^elhdi^e;2) Vl (12) 

+ ^LiLihlm + ^(LsLg + L3L2)hl7^3 + ^Li(L2/l2 V^^-s + Lshli/j^,^) 
+ '^{L2tpu,2 + L3lpu,i)hlLi . 

As shown in Section HI the vacuum structure is 

(V'e;2) = -« /Oe {'^e-l) = We , ('(/'i^;2) = Pu {lpu;l) = W^ , (r/i) = Wi and (%) = ^3 (13) 

with pe^u = il- The different choices of pe^u correspond to the different possible values of 
the subgroup index ki and k^, respectively. Since they are not uniquely fixed by the flavon 
superpotential, we keep pe^u = ±1 as parameters and check that all possibilities lead to 6*13 = 
and ^23 = vr/4. We arrive at fermion mass matrices of the form (Mi given in left-right convention) 



,,2 



Mi=\ ipeVlwi/A ipeVl -^Vd and M^ = y>^ y^w^^ -^ . 

(14) 
For 

We, Wu, wi, W3 ss eA (15) 

with e w 0.04, muon and tau lepton masses read 

m^ = V2\y2\wiWeVd/A^ ^ e^Vd and m^ = V2\yl\weVd/ J^ ^ evd , (16) 

while the electron remains massless at this stage and acquires a mass from operators with three 
flavon insertions, see Eq. ()34p . The charged lepton mass matrix can be diagonalized through the 
usual biunitary transformation (Ui and VJ) so that 

UJMiVi = diag(me, m^, m^) (17) 

with Ui given by 

/I \ 

Ui=\ ipe/V2 -ipe/V2 (18) 

V 1/^/2 1/^/2 / 

and V/ is a diagonal matrix, ensuring that the entries of Mi in the charged lepton mass basis are 
real and positive. The unitary matrix U^ diagonalizing MyMl can be written in the form 

cos (j)u sin (py e"*'''" 

Uu=\ -p^sin,/.^e^>/^/2 p^coscfiu/V^ -pu/V2 I (19) 



sin( 



^T''/^/2 cos 0^/^/2 1/^/2 



=,«7i 



^We do not list the operator Lie2ftdi/'e/A^, because plugging in the vacuum given in Ea. (|13|) renders this 
possible contribution to the charged lepton mass matrix zero. However, we discuss the operator in Section 5, 
since it leads to a non-vanishing contribution, once the shifts in the flavon VEVs are taken into account. 



with 4>u and 7jy being real functions of the entries of My. For the absolute values of the lepton 
mixing matrix we then find 

(I cos 4)y I I sin (j)i,\ \ 

|sin</.,|/^/2 I COS./., 1/^/2 1/^/2 (20) 

|sin</.,|/^/2 \cos(Py\/V2 l/\/2 / 

showing that ^23 = 7r/4 and ^13 = hold, while the solar mixing angle 612 is not fixed in our 
model and related to the angle (j)y. As one can see the form of |?7A/iV5'| and thus the result for 
the mixing angles do not depend on a particular choice of pe and p,. 

In the charged lepton mass basis the neutrino mass matrix M^ = UJMyUi reveals a texture 
zero in its (23) entry. According to the findings in [19] a consequence of the presence of this 
texture zero together with ^13 = is that the neutrino mass spectrum has to have inverted 
hierarchy and furthermore the lightest neutrino mass 7713 cannot vanish. This is in our model 
not a generic result attributed to the flavor symmetry D^ x Z5 and its breaking pattern, but is 
due to the fact that a flavon coupling to the neutrino sector and transforming as (X3, <^^) under 
(1)4,^5) is absent. A numerical analysis similar to the one found in [3] can also be performed for 
this model. The results of such an analysis are quantitatively very similar to those found in [5] 
for the case of a non-zero (23) entry of the neutrino mass matrix M'^ in a spontaneously CP 
violating framework. These include that the measure of neutrinoless double beta decay |mee| 
fulfills |?n-ee| ~ "1-3, a weak correlation between tan0i2 and m^ and rather strong correlations of 
the two Majorana phases and the mass m^ as well as a lower bound on 777,3 around 0.015 eV. 
Plots showing these results can be found in j3] in Figures 1 to 4. All results are found to be 
independent of the different choices of p^ and pe- Noting that the third neutrino mass 7773 is 
given by \y2\w'iv'^/ h? , we find that in order to correctly reproduce the light neutrino mass scale 
of around 0.1 eV a value of the cutoff scale A « 4 x 10^^ GeV for Vu ~ 100 GeV is necessary. 

Comparing our results to those found in [3] we see that the difference of the two models lies 
in the fact that in [1] the subgroup in the charged lepton sector is D2, while it is only a Z2 
subgroup in our model. As a consequence, the model in [3] leads to /i — r symmetric mixing 
through preserving certain D^ subgroups only, whereas in the model here the absence of a flavon 
transforming as X4 under D4, from the charged lepton sector is relevant for achieving ^23 = ■?r/4 
and ^13 = 0, see for details Appendix B in [3]. 

4 Flavon Superpotential 

In the following we align the VEVs of the fields V'e and ij^y correctly and fix all VEVs through 
mass parameters of the superpotential, thus avoiding flat directions in the potential associated to 
VEVs which remain undetermined. We consider the F-terms of a new set of fields, the driving 
fields, as source of the alignment conditions. These fields are like the flavons gauge singlets 
and transform in general under D/^^ x Z5 . We assume the existence of a continuous i?-symmetry 
U{1)r under which all matter superfields have charge +1, flavons and the superfields h^^d charge 
and driving fields charge +2, so that the superpotential, responsible for the alignment of the 
flavon VEVs, will be linear in the driving fields. 
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Table 2: Driving fields of the model necessary to align the VEVs of ^e and ipv and to 
fix the flavon VEVs through mass parameters of the superpotential. All these fields 
have charge +2 under U{1)r. 



4.1 Renormalizable Level 

In order to align the vacuum of ip^ we introduce one driving field, cr^, transforming as I4 under 
i^4, similar to [4], and as w under Z5. The alignment of the vacuum of V'e is achieved in a similar 
way through coupling it to a field Cg transforming as 1^3 and being invariant under Z5. We add 
a field Xe which is neutral under D4 x Z5 to allow for a coupling of mass dimension two which 
fixes the size of the VEV of V'e- The VEVs of ipi, and of the Proggatt-Nielsen fields r/i and 773 are 
deduced from the F-terms of three further fields, xl^i V^ ^'^'^ Va transforming as shown in Table 
[21 Note that also the field 7]1 allows for a coupling with positive mass dimension in the flavon 
superpotential. Apart from a'^ and cr^, responsible for the alignment of {tpe;i,2) and {'ipu;i,2), all 
other driving fields transform trivially under D4. The flavon superpotential Wfi is given, at the 
renormalizable level, asn 



W: 



}i 



«e ^e (V-e;! + V'e^) + ^l xl + be XeV'e;lV'e;2 

+ au ol [^l.^ - 11)1.2) + bu Xu ^v;i ^v;2 + Cy xlmm 



(21) 



.0^2 



+ drii (V'e;lV'^;2 + V'e;2^i.;l) + / ??3^f + M2 r?^?7l +5^/4^3 



with Oe, he, a^, b^, c^, d, f and g being complex numbers with absolute values of order one. 
Assuming that the flavons acquire their VEVs in the SUSY limit, i.e. (soft) SUSY breaking 
effects can be safely neglected, we can use the i^-terms of the driving fields to determine the 
vacuum structure of the flavons. The first two equations read 



dW 



fi 



da^ 



{<: 



We find as solutions 



+ ^e;2) = and 



dW 



fi 



dal 



(€;1 - <2) = 



(22) 



(V'e;2) = -iPe{ll^e;l) = We and (V'i.;2) = Pu{lpu;l) = W^ 



(23) 



with /)e = ±1 and py = ±1. We see that we preserve a Z2 subgroup of D4 generated by either 
BA {ki = 1) or BA {ki = 3) in the charged lepton sector, whereas in the neutrino sector the 
subgroup is either generated by B {k^ = 0) or by B A {k^ = 2). From 



dW, 



fi 



dX 







M/ + 6e'0e;lV'e;2 = O 



(24) 



^We can safely neglect the term Xehuhd, because the VEVs of all flavons are much larger than the electroweak 
scale. 



it follows that the VEV of ipe is fixed through the mass scale Mi to be 



The last three i^'-term equations, 



wl = iPe-^- (25) 



K i'u;itpu;2 + Cy rjiTj^ = , (26a) 



dWfj 

-^ = d {^e;l^u;2 + ^e;2i'u;l) + f ^l = , (26b) 



drpA 



M2rii+gr]i = {), (26c) 



allow us to determine the VEVs of V'l^j ??i and 173 as functions of the two mass parameters Mi 
and M2. For Mi and M2 being of the order eA we find that all flavon VEVs are of that order as 
well, (^e) 7^ is a necessary consequence of Eq. ()24p showing that D^ is always spontaneously 
broken. The other flavon VEVs could in principle vanish. However, if we assume that the VEV 
of one of the fields tp^^ Vi ^^^d 773 does not vanish then the non- vanishing of the other two ones 
follows. Due to the fact that all terms in Wfi are by construction linear in the driving fields the 
i^'-terms of the flavons vanish in any case for vanishing VEVs of all driving fields. Then also the 
allowed term Xe^uhd does not give rise to a //-term. 

4.2 Corrections from Non-Renormalizable Terms 

Including non-renormalizable terms with three flavons leads to corrections of the alignment 
achieved at LO. We find nine terms to contribute at this level to the flavon superpotential 
(again all possible terms involving the superfields h^ and h^ are neglected) 

AWfl = Xi a% [i^l.i + i^l.^) /A + X2 a% {lpe;l^u;l + V'e;2V'^;2) /A (27) 

+ 3^3 X^li^u-^lpu-a/^ + X4 xlm ifpe;l1pu;2 + tpe;2lpu;l) /A 

+ X5 Xef/?%/A + X6 X°r/|/A 

+ Xj r]^mi'u;l1pu;2/^ + 3^8 IJ^m'n'i/^ + ^^9 ??4^1^e;lV'e;2/A . 

Note that there is no correction at this level involving the driving field a^. We can parametrize 
the VEVs of ipe and ipu as 

(^e> = ( '^^^^^ V"""-''^ ) and m = ( ^'^(^^ +/^'^^i) 1 (28) 

\ We + dWe-2 J \ W^ + dWv-2 J 

and find for the shifts in linear expansion that 

, ~ e, (29) 

We Wy 

5we-i 7^ 6we-2 and Sw,^-i = 5wy^2 ■ 

Since the shifts of (V'z^jj) are the same, the vacuum alignment is preserved up to this level and we 
can absorb the shifts 6wy^i into a redefinition of w^. The equality of dwu-^i is due to the fact that 



at the first non-renormalizable level the superpotential AWfi does not contain terms involving 
(t[^. The shifts of the VEVs of the singlets r/i and % are of the same order of magnitude as Swe-^i 
and 5wiy-i. However, we do not mention them explicitly, because their effect on lepton masses 
and mixings can always be absorbed into a redefinition of Yukawa couplings or VEVs wi and 

5 Lepton Masses and Mixings at NLO 

In general, NLO corrections arise from two sources: (i) shifts in the flavon VEVs and (ii) oper- 
ators with two and more flavon insertions, evaluated by plugging in the LO form of the VEVs, 
which have not been considered in Section [3j As we show in the following, all such additional 
contributions change the LO results only slightly. However, their discussion is relevant, because 
such terms generate the electron mass and govern the deviations from the LO results ^23 = 7i"/4 
and ^13 = 0. We discuss all such terms which give rise to contributions up to and including e^ 
(in units of Vd and v'^/A, respectively) as well as for the (13) element of the charged lepton mass 
matrix also corrections of order e^v^. 

Considering the charged lepton sector we find the following additional operators involving 
two and more flavons which contribute to the (23) and (33) elements of M; (order one coefficients 
are omitted and operators are not written in D4 components) 

1 1 o 

j2 ^Deshdip^rjs + ^ Loe^hdipe (30) 

where the last term actually gives rise to two independent contributions. In a similar way the 
(22) and (32) elements get corrected through 

-p Loelhd'ilft + -j^ Loe'ihdiperil + j^ LDe2hdil^uViV3 (31) 

with the first term being responsible for two independent contributions. At the same time these 
entries also receive contributions from the shifts in the VEV of V'e which we indicate by □ 

— LoelhdS^e + -r2 LDelhddipem ■ (32) 

The (21) and (31) elements are generated only through three fiavon insertions of the form 

-3 Loelkdipei'l + -Tg LDelkd-ipv^l + -r^ LDelhdipemVs ■ (33) 

Note that the first operator gives rise to three independent terms. Similarly, the (11) element 
arises from 

-r^ Lielkdipe-^vm . (34) 

The generation of the (12) and (13) elements is somewhat special, because the lowest order 
operators which could give rise to non-zero (12) and (13) elements are 

-^ Lie^hdipl and -^ Lieg/i^V^ . (35) 



^As mentioned in Section [4. 21 the shift in the VEV of the field ri\ can be absorbed into the Yukawa couphngs 
or the LO VEV itself and thus its contribution is not displayed. The same holds for 773. 



However, plugging in the LO result for the VEVs of ipe and ipi, we see that these operators give 
no contribution. Taking into account the shifts arising at a relative order e we find that then the 
(12) element is generated, whereas the (13) element still vanishes. Operators with three fiavons 
also contribute to the (12) element at this level 



The (13) element only originates from operators involving four fiavons 



(36) 



(37) 



Again, the last operator allows for two independent contractions. At this level we expect that 
the (13) element also receives a contribution from subleading shifts in the VEVs of ipu which 
we however do not calculate, because the specific form of the (13) element is not relevant for 
the analysis of lepton masses and mixings. Thus, we can parametrize the charged lepton mass 
matrix including NLO corrections as 



/3|62 



Ph^ 



Pie' 
Mi=\ /3|e2 ip,ale + Ple' ipeaf+^^e \ ev^ 
/3|e2 a|e -af 



(38) 



taking into account the sizes of fiavon VEVs and of their shifts as given in Eq. (|15p and Eq. (|29p . 
Through re-phasing of right-handed fields we can make af, a\ and /3| real and positive. The 
other parameters, apart from e, are in general complex numbers with absolute values of order 
one. Note further that the parameters a\ 2 are up to corrections of order e determined by the 
LO operators given in Eq. (fT2|) . 

We find as result for the charged lepton masses 



me. 



l\e- + 0(e^)) Vd, m^= (V2ale^ + 0{t^)) va and m^ = (V2ale + 0{e 



Vd 



(39) 

and thus confirm earlier statements about the size and origin of the electron mass in our model. 
The unitary transformation applied to the left-handed charged leptons in order to diagonalize 
Ml is, up to the first correction in e in each matrix element, of the form 



I i-OM)\^ 



Ui 



2a% 



^1 



-ip^ 



2a% 



\ 



2a% 



1 

1 

V2 



ipe + 

l + ip, 



^^gS 



:2af ' 



ipe + 



'71 \''^'^ ^ ^^ 



\ 



J 



(40) 



We analyze the corrections to the neutrino mass matrix in a similar way and find that the 
following operators contribute to the (11) element 



^ LiLihlrJl + ^ LiLihl'4)lr]i 



and that corrections to the (23) element come from 

LoLDhlni + —J {LDLD)ihl{ij'^ri3)i + — (LDLD)i/i^(^eV'i^??i) 



A4 



A4 



A4 



(41) 



(42) 
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Here we denote by (•••)! the contraction to a trivial singlet 1^ of D4. All these contributions 
can be absorbed into the LO result. The relation between the (12) and (13) element is disturbed 
at a relative order e^ through the terms 

-^ LiLohlil^lil^y + ^ LiLohli^eVi ■ (43) 

Note that the first operator generates three independent contributions. Since the shift of the 
VEVs of Tpi, is at lowest order aligned with the LO result, see Eq. (p9|) . we do not encounter 
any correction to the LO relation between the (12) and (13) elements at the relative level of e. 
Among the corrections of order e'^ (in units of v'^/A) we also expect corrections due to possible 
deviations from {ipy^i) = {il^u;2) at the level of e^A. The (22) and (33) elements of Mj, receive 
two types of corrections. The first one still preserves the n — t symmetric structure of the LO 
result and generates entries of order e^f^/A 

^ {LDLD)3hl{tplm)3 + ji {LDLDhhl{iJei^^rii)3 ■ (44) 

The second type breaks jU — r symmetry and stems from very similar operators 

-^ {LDLD)4hl{'ipl'n3)4 + ^ {LDLD)Ahl{ll)ei^vm)i ■ (45) 

As one can see the only difference is the contraction to another D^ singlet; in the first case we 
contract to a 1^3, indicated through (• • • )3, whereas in the second case the contraction is to I4) 
denoted by (• • • )4. Also these fj, — t symmetry breaking contributions arise at the level e^f^/A. 
Thus, at NLO M^ takes the form 



M,= I ag+/3re' (3!^e^ a^ | e ^^ . (46) 






A 



Note that the parameters 0123 and /5i2 3 ^^^ complex numbers with absolute values of order 
one. Similar to the charged lepton sector, the parameters Oi 2 3 coincide at LO in the expansion 
in e with those given in Eq. (jl2p and Eq. (jl4p . 

Recomputing the mass spectrum of the neutrinos we find that all masses are corrected by 
terms of relative order e^. Relevant for calculating the deviations from maximal atmospheric 
mixing and ^13 = is the form of the eigenvector associated to the third neutrino mass 777,3. We 
find that its LO form, see Eq. ()19p . receives corrections only at order e^. Thus, the results for 
the mixing angles ^13 and ^23 at NLO are dominated by corrections coming from the charged 
lepton sector, see Eq. iflO]) . and we get 

|C/e3|=sin%«^e, (47) 

The solar mixing angle 9i2 whose exact value is not predicted in our model, but only given in 
terms of the parameters of the neutrino mass matrix at LO, undergoes also corrections of order 
e. However, these are not of particular interest. 
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6 Conclusions 

We have presented a SUSY D4 model which leads to maximal atmospheric mixing and vanishing 
^13, while keeping the solar mixing angle 612 undetermined, but expected to be in general 
large. These predictions originate from the non-trivial breaking of D4 to distinct Z2 subgroups 
generated by BA , k = 0,1,2,3, in the charged lepton and neutrino sectors. As we have 
shown, it does not matter to which of these Z2 subgroups D4 is actually broken. It is enough 
to achieve that the index ki in the charged lepton sector is odd, whereas k,^, the index in the 
neutrino sector, is even. At the same time the hierarchy among charged leptons can be naturally 
produced. Apart from D4, responsible for the ^i — t symmetric mixing pattern, we employ a 
cyclic symmetry Z5 in order to separate the flavons coupling to neutrino and charged lepton 
sectors at LO as well as to achieve appropriately the mass hierarchy me <C rn,^ <^ m-r- Due 
to a restricted choice of flavon fields we find that neutrinos have inverted mass hierarchy with 
the lightest neutrino mass m^, > 0.015 eV. These properties together with non-trivial relations 
among TO3, the observables |r?T-ee| aud tan0i2 as well as the Majorana phases are shared with 
a similar SUSY L>4 model, found in 0]. In the latter, however, the hierarchy among charged 
leptons could not be reproduced without fine-tuning. This is due to the fact that right-handed 
charged leptons also transform as l^-|-2 under D^. This problem is solved in the present model 
by assigning them to three (distinct) one-dimensional representations of -D4. The result oi jjl — t 
symmetric lepton mixings remains untouched. 

A particular feature of our model is that the flavon superpotential not only leads to the 
vacuum alignment which is needed in order to arrive at the predictions ^23 = '?r/4 and ^13 = 
0, but also fixes all flavon VEVs through two mass parameters. Mi and M2. In this way, 
it is natural to have the same order of magnitude for all flavon VEVs and flat directions in 
the flavon (super)potential, associated to free parameters among the VEVs, are avoided. In 
contrast, in (almost) all models of this type these are generically present. The existence of 
terms with couplings of positive mass dimension also ensures the spontaneous breaking of the 
flavor symmetry D4 and so destabilizes the trivial vacuum in which all flavon VEVs vanish. 

We have carefully studied NLO corrections arising from higher-dimensional operators involv- 
ing several flavons in the Yukawa sector as well as in the flavon superpotential. The latter induce 
in general shifts in the flavon VEVs. In the particular case we discussed here these operators 
do not disturb at NLO level the LO vacuum structure of the flavons coupling dominantly to 
the neutrino sector. We eventually find that all corrections in the neutrino sector are of relative 
order e^. In contrast to this, the results in the charged lepton sector get corrected at a relative 
order e. For this reason, the deviations from ^23 = 7r/4 and ^13 = of order e stem from the 
charged lepton sector only. This is a further feature which distinguishes the present model from 
the one of |3] , in which deviations from ^ — t symmetric lepton mixing are governed by correc- 
tions from the neutrino sector and which at NLO lead to ^23 — 7r/4 being much smaller than ^13. 
In the present model NLO corrections are furthermore relevant to generate the electron mass 
of order e^Vd- As frequently happens in such models, the mass of the tau lepton is generated 
through a non-renormalizable operator so that small and moderate values of tan /? are preferred. 

Several models leading to ^23 = 7r/4 and ^13 = have been considered before in the literature. 
The model closest related to the one illustrated is discussed in [4j. However, as mentioned, this 
model needs some fine-tuning to get the correct charged lepton mass hierarchy. The model [6], 
on which |3] is actually based, is in contrast non-SUSY and thus incorporates no solution to 
the hierarchy problem. Nevertheless, in this model it has been shown [6l[20] that m,^ <C rrir 
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can arise from a softly broken (additional) symmetry, whereas nie being small can be attributed 
to a small VEV of one of the Higgs doublets present in the model. A model constructed in 
the same spirit as [6] can be found in [7] having as flavor group ^3 x Z2 (recall that S3 is 
isomorphic to -D3). Its results are very similar to those of the model in [6j, apart from fewer 
constraints on the neutrino mass spectrum. A neat example of a minimalistic SUSY model 
leading to fi — r symmetric mixing is given in [T^- The flavor group is 53 x Z3. |£| The structure 
of the model is very similar to ours and in the neutrino sector also a Z2 subgroup of the 
same type as here is preserved by appropriate flavon VEVs. Yet, in the charged lepton sector 
an alignment is employed which completely breaks 53, but efficiently generates the charged 
fermion mass hierarchy. The alignment is studied in detail in |12) . but, as usual, two parameters 
among the flavon VEVs remain undetermined, thus giving rise to flat directions in the flavon 
(super)potential. 

In summary, we have constructed a SUSY D4 model, which predicts /i — t symmetric lepton 
mixing through breaking Z)4 to distinct Z2 subgroups in the charged lepton and neutrino sectors, 
respectively, and at the same time accommodates naturally the mass hierarchy rUe <C m.^ <C rrir- 
Furthermore, all flavon VEVs are fixed through mass parameters of the superpotential. Thus, 
the problem of free parameters among the flavon VEVs is avoided, which is often met in models 
of such type. On the basis of this model it might be very interesting to consider an extension to 
the quark sector. As has been shown in [^, also the Cabibbo angle might arise from a non-trivial 
breaking of a dihedral group to distinct Z2 subgroups in up and down quark sectors and so its 
size might have a similar origin as the prediction of 623 = vr/4 and ^13 = in the lepton sector. 
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